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ABSTRACT 



Spiral arms, if they are massive, exert gravitational torques that transport 
angular momentum radially within galactic disks. These torques depend not on 
the pattern speed or permanence of the arms but only on the nonaxisymmetric 
mass distribution. Hence the torques can be measured from photometry. We 
demonstrate this using gri CCD data for MlOO (NGC 4321). Since we find 
consistency among the three bands, we believe that dust and young stars in the 
arms do not seriously bias our results. If the present epoch is representative, the 
timescale for redistribution of angular momentum in MlOO is 5 — 10 Gyr, the 
main uncertainty being the mass-to-light ratio of the disk. 
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1. Introduction 

It has long been recognized that spiral structure may be associated with secular 
dynamical changes in galactic disks. The anti-spiral theorem of Lynden-Bell & Ostriker 
(1967) states that strictly steady state spiral modes of linear amplitude cannot exist 
without dissipation, except where they are superpositions of degenerate nonspiral modes. 
Nonlinear selfconsistent models of spirals have been sought [Contopoulos & Grosb0l (1988), 
Contopoulos (1990) and references therein]. Lynden-Bell & Kalnajs (1972) have shown, 
however, that massive spiral arms transmit gravitational torques. If the arms trail, then 
the torques transport angular momentum outward through the disk. Linear or nonlinear, 
dissipative or conservative, a rigorously steady spiral pattern therefore requires sources and 
sinks of angular momentum. But is this a problem only of principle; in other words, how 
significant are spiral torques compared to the ages and angular momenta of their galactic 
hosts? 

On the basis of extensive photographic photometry, Schweizer (1976) argued that red 
(or rather orange) light traces the old stellar disk in external spirals. If this is true then 
spiral arms represent an enhancement in the surface density of mass. Since Schweizer's 
influential study, many workers have Fourier analyzed photometry of spirals in red or 
near-infrared light to measure spiral amplitudes: for example, lye et al. (1982); Considere & 
Athanassoula (1988); and Elmegreen et al. (1992). Although these and other studies made 
quantitative measurements of the strength, coherence, and multiplicity of the optical arms, 
the concern has persisted that the relationship between light and mass may be distorted 
by dust and by bright young stars formed in the arms. Recently, the advent of solid-state 
infrared detector arrays has allowed photometrically precise images at J, H, and K, where 
dust (probably) and young stars (arguably) are much less problematic than they are at 
visible wavelengths. In a detailed study of M51, Rix (1993) and Rix & Rieke (1993) have 
made a good case that K light traces mass — or to be more cautious, that the extinction at 
if is ^ 10% and that the nonaxisymmetric structure at K appears not to be dominated by 
stars whose age is less than the interarm crossing time. These authors also found that the 
/-band surface brightness was very faithful to that at K . Quillen et al. (1994a) and Quillen 
et al. (1994b) have used JHK photometry to reconstruct the nonaxisymmetric potentials of 
the stellar bars in NGC 4314 and NGC 7479. 

The object of the present study is to measure spiral-arm torques directly from 
photometry. This seems not to have been done before. The torques depend only upon the 
mass distribution, and not upon the persistence or pattern speed of the arms. One needs a 
photometric tracer of the mass. Although a dark-matter halo that is not at all well traced 
by the light is probably present and may dominate the rotation curve, this is not a serious 
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concern for the torque measurement. Dynamical stability requires that the halo be much 
hotter kinematically and thicker geometrically than the visible stars [Ostriker & Peebles 
(1973)], so the halo probably does not participate in spiral structure. 

To translate the torque into a secular evolution rate, one must assume that the torque 
measured today is typical of the galaxy at other times, although one need not assume that 
the spiral pattern is steady. Also, the torque has to be compared to the angular momentum 
of the stars in order to estimate how rapidly their orbits change, which requires knowledge 
of the rotation curve. 

To illustrate the feasibility of the torque measurements, we have applied our methods 
to the bright and well-studied Virgo spiral MIOO. The data are described in §2, together 
with our assumptions about dust obscuration, mass-to-light ratios, and disk thickness. The 
method of analysis is laid out in §3, together with tests on artificial data. Results for MIOO 
are given in §4, and §5 sums up. 

2. MIOO 
2.1. The Images 

Images of MIOO were taken from A Catalog of Digital Images of 113 Nearby Galaxies 
[Frei et al. (1995), henceforth FGGT]. Observations were made on the 1.5 meter telescope 
(P60) of the Palomar Observatory on the night of 4 May 1991 with the Wide Field Prime 
Focus Universal Extragalactic Instrument. Images were taken in the g, r, and i bands of 
the Thuan-Gunn photometric system [Thuan & Gunn (1976), Wade et al. (1979)]. These 
filters are centered at 0.50 /x, 0.65 /x, and 0.82 /x, respectively. See Frei & Gunn (1994) and 
FGGT for details. 

The field of view is 16' X 16', and it is projected on an 800 X 800 Texas Instruments 
CCD chip in the Cassegrain focus, yielding a scale of l."19 pixel"-^. Using Hubble Space 
Telescope observations of cepheids, Freedman et al. (1994) have estimated the distance of 
MIOO as 17.1 lb 1.8 Mpc. Thus the projected size of a pixel is 99 it 10 pc. 

A uniform sky level was subtracted from the images and foreground stars were removed 
[see Frei (1995)]. The photometric calibration by FGGT used observations of seven standard 
stars. We have converted the calibrated data in gri to the corresponding Johnson bands 
VRI according to the transformations of Frei & Gunn (1994). After transforming the 
colors and rectifying the images (see below), we have compared the azimuthally-averaged 
surface brightness profiles and total magnitudes with the photometry of de Jong & van der 
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Kruit (1994). Thus, whereas we find mj = 8.38, tur = 8.99, and my = 9.47 (errors not 
estimated), de Jong & van der Kruit (1994) quote m/ = 8.56 it 0.11, ruR = 9.01 it 0.10, 
and mv = 9.56 it 0.08. All magnitudes refer to the light within R25, the radius of the 
25th-magnitude isophote (at B). From the data of FGGT, we estimate R25 = 221" or 
21.9 kpc, in excellent agreement with the value 222" cited by the RC3 [de Vaucouleurs et 
al. (1991)]. Therefore, we believe that the photometry and color transformations available 
to us are sufficiently accurate for our purposes. 

The images were rectified to correct for inclination. This is important, because tests 
on artificial data indicate that the torque estimate is sensitive to rectification errors (§3). 
There is no doubt that MlOO is seen rather close to face-on, but published estimates of the 
position angle and inclination vary considerably. Having made the attempt ourselves, we 
are skeptical of estimates based on photometry alone. The ellipticity caused by inclination 
is second-order in the inclination [i) when i is small; and it is entangled with the intrinsic 
m = 2 spiral structure. The observed velocities are first-order (oc sin 2); and although spiral 
structure may disturb the velocity field, the perturbations are also proportional to sin 2. It 
is relatively easy to identify the kinematic axes. Therefore we have adopted the position 
angle 158° determined by Arsenault et al. (1988) from Ha Fabry-Perot data, which is in 
good agreement with other kinematic measurements [Rubin et al. (1980): 155°; Warmels 
(1988), Guhathakurta et al. (1988): 153°]. Some photometric determinations are almost 
orthogonal, however: de Vaucouleurs et al. (1991): 30°; Grosb0l (1985): 58 ±4°. 

We determine the inclination by requiring that the rotation curve agree with the 
Fisher- TuUy (FT) relation. With the total magnitudes tur and mj cited above and the 
reddening corrections derived in §2.2, the relations of Pierce & TuUy (1992) predict a 
circular velocity Vc = 270 it 22 km s""^, where the error is the sum in quadrature of Freedman 
et al. (1994)'s uncertainty in the distance modulus, the 0.2 mag difference between our mj 
and that of de Jong & van der Kruit (1994), and the 0.1 mag intrinsic width of the FT 
relation in R and /. Comparing this with the HI data of Guhathakurta et al. (1988), we 
find 2 = 27° ± 2.5°. 

We have used the central 512 X 512 pixels (609" X 609") of each field for our calculations. 
Grayscale plots of the original 2, r, g images and the rectified g image are given in Fig. 1. 
The g image shows more small-scale structure than the 2, presumably as a result of dust 
and recent star formation. The straight line in Fig. Id shows the major axis (line of nodes) 
we have adopted. The small companion seen at the top of the image has proved to have 
negligible infiuence on the measured torque. Two prominent arms are clearly seen. They 
are more open in the rectified image than in the original. 
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2.2. Corrections for Galactic extinction and internal absorption 

Since we want to use light as a tracer of mass, it is important to correct the data for 
the effects of dust. 

According to Burstein & Heiles (1984), Galactic extinction of MlOO in the B band 
is only 0.04mag. For a standard extinction curve [Savage & Mathis (1979)] that 

is compatible with Burstein & Heiles (1984)'s assumption Ab/E[B — y) = 4, we find 

% 0.031, A% % 0.023, Af % 0.014 for the Galactic extinction at V, R, and /. These 
values are so small as to be essentially negligible compared to uncertainties in the stellar 
mass-to-light ratios, but we correct for them anyway. 

The internal absorption within MlOO can be estimated in several ways. TuUy & 
Fouque (1985) have derived an empirical formula for internal absorption, which is based on 
the trend in B — H color with inclination in a sample of late-type spirals. Their formula 
predicts < >^ 0.30 mag for i = 27° (and 0.27 mag for i = 0°). The reduction in 
surface brightness due to the dust, which is what these numbers measure, is not directly 
proportional to the extinction optical depth of the disk. The relationship between surface 
brightness and extinction depends on the relative thicknesses of the star and gas layers, 
and on scattering. Ignoring these fine points (because the final correction is small), we use 
the standard extinction law to derive [Ay,A]^,A^j) = (0.23,0.18,0.11) mag. We can check 
this against Freedman et al. (1994)'s value for the mean extinction of 20 Cepheids in MlOO, 
as derived by comparing the period-luminosity relations in V and /: these authors have 
obtained (Av) = 0.15 it 0.17. Finally, Xu & Buat (1994) have used the far-infrared emission 
of 135 spirals to deduce mean and median total face-on optical depths at B of 0.60 it 0.04 
and 0.49, respectively. The surface-brightness reduction corresponding to these optical 
depths depends on details of the radiative transfer but is probably about half as large; 
thus, (Ab) ~ 0.3 mag. We regard these three ways of estimating the internal absorption as 
roughly consistent. 

Our adopted values for the total corrections to the surface brightness, accounting for 
both Galactic extinction and internal absorption, are Ay = 0.26 mag, Ar = 0.20 mag, and 
Aj = 0.12 mag. Hence the redenning-corrected total magnitudes of MlOO are estimated to 
be mi = 8.26, tur = 8.79, my = 9.21. 

2.3. Stellar Mass-to-Light Ratios 

The mass-to-light ratio of the stellar population in MlOO (or any other spiral) is the 
greatest uncertainty in estimating spiral-arm torques from photometry alone. We have used 
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published population-synthesis models for this purpose. 

Following the usual convention, we define [M/L)x, the mass-to-light ratio in band X, 
as the mass of the stellar population in units of Mq divided by the number of solar-type 
stars that would produce the observed light in band X, after correcting for absorption. 
Thus if M is the mass of the population, Mx its absolute magnitude in band X, and Mx,© 
the corresponding absolute magnitude of the Sun, then 

^=(^\ 1oO.4(Mx,0-Mx) (I) 

\l)x ^ ' 

The models of Arimoto & Jablonka (1991) as interpreted by Jablonka & Arimoto 
(1992) predict [M/L)h = 2.15 for Sc galaxies, including the mass in associated gas and 
stellar remnants, but not halo dark matter. This requires assumptions about the slope 
and extent of the initial mass function and about the history of the star-formation rate. 
The mass-to-light ratios at V, R, I can be derived from [M/ L)h using the predicted colors 
{V - H)sc = 2.196, {R - H)sc = 1.567, and (/ - H)sc = 0.989, and the colors of the 
Sun [Worthey (1994)], {V - H)q = 1.47, {R - H)q = 1.11, and (/ - H)q = 0.76. Thus: 
{M/L)v = 4.20, {M/L)r = 3.28, and {M/L)i = 2.65. 

For comparison, the single-burst models of Worthey (1994) predict [M/L)i 
monotonically increasing from 1 to 3 with ages from 1.5 to 10 Gyr. According to these 
models, the dependence on metallicity is very weak in the / band. For populations of age 
5 Gyr and solar metallicity, [M/L)v =3.1 and [M/L)r = 2.6. These bands are more 
sensitive to age and considerably more sensitive to metallicity than the / band. All of 
Worthey's mass-to-light ratios are based on the initial mass in stars less massive than 2Mq. 
We refer the reader to the papers of Arimoto & Jablonka (1991) and Worthey (1994) for 
further details. 

Both the absorptions cited in the last subsection and the mass-to-light ratios cited in 
this one are mean values for the entire disk. The absorption is larger than the mean in 
the arms, however, where the gas and dust concentrate. Furthermore, the arms contain 
a higher proportion of young, luminous stars. Both of these effects, which have opposing 
infiuences on the effective mass-to-light ratio, are expected to be more severe at V and R 
than at /. When calculating the surface mass distribution from the photometry, we have 
not corrected for these variations in mass-to-light ratio. The extent to which our conclusions 
are compromised by mass-to-light variations can be gauged by comparing our results for 
the mass distribution and the torque estimated separately from the data in the three bands. 



2.4. Thickness of the Disk 
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Infrared photometry of edge-on spirals indicates that the light profile perpendicular to 
the disk is constant with radius and well fit by an exponential law, 

P(^) = '''"'^If'-' . (2) 

with a vertical scale height that is of order one tenth the horizontal disk scale length, hr. 

In their infrared study of IC2531, Wainscoat et al. (1989) found ^ hr/12. In a 
similar study of N5907, Barnaby & Thronson (1992) fit a sech^z / h^) profile — suggested by 
van der Kruit (1988) as a compromise between the exponential law (2) and the dynamically 
simpler sech.^{z / hz) profile — and found hz/hj. ^ 1/9.3. Kuijken & Gilmore (1989) fit their 
counts of K dwarfs in the solar neighborhood to a double exponential: 82% of the stars are 
in a component with hz = 249pc, and the rest have h'^ = l.OOkpc; they assume hr = 4.5kpc, 
so hz/K ^ 1/18 and h'^/hr ^ 2/9. 

For the purpose of computing the torque, we have assumed the profile (2) with 
hz = K/12. 



3. Recipe for the Torque 



In this section we explain our methods for calculating the torque from the photometric 
data. We assume for now that we know the mass-to-light ratio and can convert the surface 
brightness, o", to the stellar mass per unit area, E. 

Let (r, 6, z) be cylindrical coordinates such that r is the distance from the rotation axis 
of the galaxy, z is the height above the equatorial plane, and 6 is an azimuthal angle. The 
z component of the gravitational torque across radius R, exerted on the outer (r > R) part 
of the galaxy by the inner (r < R) part, is 

T{R) = G [ d?x,p{x,) f d'x^pix^)^-^^^^^ 

J J \Xi — X2 

T-l>ii T-2<ii 

oo 

= Gjde C{e) J d^r.^r,) J ^^^^^(^0 [|^^ L'^.y^ ^^js/^ > (3) 

— oo T'l>ii r2<R 

where Xi is the vector corresponding to coordinates [r,6,z), r is the projection of this 
vector onto the equatorial plane, p{x) is the stellar mass per unit volume, and S(r) is the 
mass per unit area. In the second line, we have assumed the disk has a constant vertical 
structure, so that p[r,z) = 'Z[r)p[z) with / dz p{z) = 1, and we have introduced the 
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vertical autocorrelation 



oo oo 

C(e) = J dz p{z)p{z + e), J de C{e) = 1. (4) 



It is convenient to replace the integration over e in Eq. (3) by evaluation at e = Coj 
for an appropriate choice of Co- Taking Co = is equivalent to treating the disk as though 
it had zero thickness. In fact the disk has some characteristic thickness ~ 2/p(0). In 
the limit kR ^ 1, T[R) can be decomposed into independent contributions from each 
Fourier component E^(r) oc exp[ik ■ r). For fixed amplitude, the contribution from 
any such component is a decreasing function of the thickness: if kh^ <C 1, then the 
thickness hardly matters, but if kh^ ^ 1, the reduction in the torque is substantial. Thus 
if p{z) has the exponential form (2) then it can be shown that the reduction factor is 
(1 + IAj/i^I)"-^. If one replaces C(e) 6[e — Co), however, then the reduction factor becomes 
exp( — |A;eo|)- Our prescription for Co is that the exponential function should match the true 
reduction factor when the latter falls to 1/2. In the case of the profile (2), this leads to 
Co = hjn{2) % 0.693K} 

We decompose the surface mass distribution into its Fourier components over azimuthal 
angle 6 at each radius 

oo 

T,{r,6) = So(r) + ^ Ecm cos(m^) + E^^ sin(m^). (5) 

m=l 

The softened newtonian potential can be expanded in a series of modified Laplace 
coefficients: 

oo 

^ Fm{ri,r2;eo)cosm{ei - 62), (6) 



irl + rj - 2rir 2 cos{9i - 62) + el m=o 
in which [Prasad (1932)] 

Fm{rur2;eo) = — ^= Qm-1/2 • (7) 

Here Qm- 1/2(2;) is a Legendre toroidal function of the second kind and Sij is the Kronecker 
delta. 



■^The isothermal profile p{z) = sedi^^z/ hz)/2hz is sometimes used. In this case the 
reduction factor is [1 — ^ + Y2C^'0(2, 1 + V2C)]) where ( = khz and 'ip{n, z) is the polygamma 
function, anci Sr. ~ 0.7812/1,. 
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We can now integrate Eq. (3) over azimuthal angles to obtain 

oo 

T{R) = Y: ^m{R), (8) 

m=l 

poo pR 

Tm{R) = GiT^m ridri c?r2 Fm(ri, rz; £„) [Ssm(n)Scm(r'2) - Ssm(r'2)Scm(n)] • 

JR Jo 



The torque vanishes if the surface density is bilaterally symmetric around any horizontal 
axis, since if we were then to orient our coordinates so that the symmetry axis lay along 
^ = 0, then Esm(r) = for all m and r. Whatever the shape of the galaxy, there is 
no contribution to the torque from any m = component of the density, such as an 
axisymmetric dark halo. 

The torque due to a spiral arm of the form 

T,m{r,e) = T,mo{r)cos{me + (j)m{r)) (9) 

approaches asymptotically 



(/32 + m2)3/2 ^' ^ din', 

in the limit 



(10) 



r=R 



m.^fi' + ,n'»ra^[^l.^,^j. (11) 

Constant /3 gives a logarithmic spiral of pitch angle tan~-^(m//3). At a given k, the torque 
is maximized when /3 = m, which corresponds to a very loosely-wrapped spiral with pitch 
angle = 45° . 

The contribution of the disk to the rotation curve can be expressed in terms of the 
m = component of the surface density: 

V.'4^sk{r) = -2^Gr— J ^o{r')Fo{ry-eoydr'. (12) 



The actual rotation curve, V^, includes contributions from the bulge, spheroid, and dark 
halo, and corrections for asymmetric drift. 

Given Vc{r), we estimate the angular momentum of the disk interior to radius R: 

l-R 

J{R) = 2'K rVc{r)Y.o{r)rdr. (13) 

^0 
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The angular momentum within R changes because of the torque across R and because of 
changes in the mass within R: 



-iR,t) = -T(R,t) + ——(R,i). 
Hence, a characteristic evolution time at radius R is 



(14) 



tev{R) 



J{R,t) 



(15) 



We have tested the above procedures against artificial data for m = 2 logarithmic 
spirals with two radial profiles: a power-law disk [S2('") oc r~^] and an exponential disk 
[E2(r) oc exp(— r//ij.)]. Like the real images, our artificial data is binned into 512^ pixels. 
The asymptotic expression Eq. (10) agrees very well with the integral formula Eq. (8) 
for the exponential density profile; the relative error is a few percent when the "image" 
subsumes about four scale lengths (as do the real data — see §4). But the power-law disk 
is very sensitive to truncation of the image. Whereas the correct value of T[R) should 
increase linearly with radius, assuming unlimited radial extent and zero thickness, the 
torque calculated from the artificial data oscillates with increasing amplitude around the 
correct value. 

Substantial errors occur in the torque calculated from Eq. (8) if the image is incorrectly 
centered or rectified. An incorrect center gives rise to a very strong m = 1 harmonic, and 
the torque oscillates strongly with radius around its correct value. Improper inclination 
amplifies the torque vastly and also produces some variations with radius. This does not 
seem to be a severe problem in the case of MlOO, as we discuss in §4.3. 



4. Results 



4.1. Photometric results 



We have applied the techniques described in §3 to MlOO. We decompose the surface 
light distribution into its Fourier harmonics, o"cm('") and o"sm(r), for m < 10. Since the 
Fourier components (see Fig. 4 below) decrease rapidly with m, and their contributions 
to the torque decrease even faster [cf. Eq. (10)], the first ten harmonics are more than 
sufficient to reproduce the relevant features of the surface density. Fig. 2 shows images 
reconstructed from these Fourier components in i and g. The reconstructions are faithful to 
the originals, but smoother. 
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The axisymmetric surface-brightness profile, cro(r), is shown in Fig. 3 in semilogarithmic 
form. In this and all subsequent figures, we have transformed the photometric data into 
the Johnson VRI bands, as described in §2.1. A compact component, presumably the 
bulge dominates the light within 20". Between 60" and 160", the disk can be approximated 
by an exponential with scale length hr ~ 64" (5.35 kpc). The three bands are parallel in 
this diagram to better than 0.3 mag between 10" and 190"; at larger radii / declines more 
rapidly than V or R. 

Profiles of the nonaxisymmetric Fourier components are shown in Fig. 4. The moduli 
and phases are very consistent across the three photometric bands within 200", especially 
for m > 2. 

Figure 5 shows the azimuthal variation of surface brightness at five selected radii. 
Following Rix (1993, Fig. 8), we have subtracted all of the odd Fourier harmonics. While 
it is clear that these profiles are not just m = 2 sinusoids, the differences among the three 
bands are rather small. The maximum brightness contrast between arm and interarm 
regions is about 1.2 mag, or a factor of 3, in the outer parts of MlOO; thus we are in 
qualitative agreement with Elmegreen et al. (1989), who found contrasts as large as a factor 
of 4. 

4.2. The rotation curve 

We have used Eq. (12) to calculate the contribution of the disk to the rotation curve; 
the results for all three photometric bands are shown in Fig. 6. A constant vertical scale 
height hz = 446 pc has been assumed even in the "bulge" (cf. §2.4). Because of the close 
similarity of the axisymmetric light profiles (Fig. 3), the three rotation curves are similar in 
shape, but with Jablonka & Arimoto (1992)'s mass-to-light ratios, they differ in amplitude. 
Therefore, we have fixed [M/ L) in each band so as to minimize the differences among the 
three curves in the range 60" < R < 200" region) and to cause them to peak at about 
270 km s""^ (§2). The required ratios of dynamical mass to light are 2.85, 3.49, and 3.80 in 
/, R, and V , respectively (in solar units, after correction for dust). 

The dynamical mass includes the mass in the dark halo within the radius where 
the rotation curve peaks (-Rmax ~ 160"). Recall that the corresponding stellar (M/iv)'s 
predicted by the population-synthesis models are 2.65, 3.28, and 4.20 (§2.3). Thus the disk 
would appear to account for a fraction / ^ 0.93 (based on /) of the mass contributing to 
the observed part of the rotation curve. Although the predicted colors are in fair agreement 
with the observations — e.g. [V — I)sc = 1-21 versus [V — /)mioo = 0.95 (or 1.09 before the 
reddening correction) — -the models of Arimoto & Jablonka (1991) and Worthey (1994) may 
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overestimate the mass in stellar remnants and very-low-mass stars, which would contribute 
little light in any of our bands. Rotation-curve fitting [Persic & Salucci (1990)] indicates 
/ ^ 0.5, a value typical of bright spirals. The bar-instability criterion of Ostriker & Peebles 
(1973), if applicable, also requires / ^ 0.5. 

4.3. The Torque and the Evolution Time 

The torque computed from all three bands is shown in Fig. 7. Light has been 
translated to mass according to the dynamical ratios above, reduced by a factor / = 0.93. 
Because it is quadratic in the surface density, the torque scales as The agreement 
among the three bands is clearly very good. As shown by the lower panel, almost all of the 
torque is created by the m = 2 Fourier components. The torque peaks at -Rpeak ~ 130", 
or about 2 scale lengths. It falls off at larger radii because the disk as a whole quickly 
fades, although the arm-interarm contrast remains strong. The local rate of accumulation 
of angular momentum is proportional to —dV/dR, so the disk interior to -Rpeak is (at 
present) losing angular momentum, while that exterior is gaining it. The slightly negative 
values of T at R ^ 80" may be due to errors in centering or rectifying the image. It 
is also possible that we are picking up a weak leading wave (assuming that the spiral 
structure is predominantly trailing), whose presence would be required in steady-state by 
the swing- amplifier mechanism [e.g. Toomre (1981)]. 

As we indicated in §3, the measured torque curves are sensitive to the rectification of 
the galaxy. The upper panel of Fig. 8 shows the results of tests on an artificial galaxy with 
an exponential profile and an m = 2 logarithmic spiral. The hatched region shows the range 
swept out by the "measured" torque when we varied the assumed position angle by ±10° 
and the inclination by ±5° around the true values. The data for MlOO are less sensitive 
to the rectification parameters. In the lower panel of Fig. 8 the /-band torque is plotted 
for our preferred values, i = 27° and PA = 158°. The hatched region is the range swept 
out when the PA is varied by ±10° and i by ±5°. Our tests on artificial data show that 
the sensitivity to rectification increases as the spiral arms become more tightly wrapped. 
MlOO is a relatively loose spiral, hence relatively insensitive. The artificial spiral used in 
constructing the upper panel was chosen to have a pitch angle comparable to that of MlOO 
(/3P.1.7). 

Figure 9 depicts the evolution time Eq. (15) for the three bands. Because of the way 
tev(-R) has been defined, it must diverge at large radii, where the torque vanishes, and it 
isn't reliable at small radii, where the spiral structure is poorly defined. Therefore, the most 
interesting features of these curves are their minima, which occur at about two disk scale 
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lengths: tev.min ~ 6 Gyr. In terms of the mass fraction / introduced above, J{R) oc / for 
a fixed amplitude of the rotation curve, so tgv oc f~^. Thus if / ~ 0.5 instead of 0.93, the 
values shown in Fig. 9 should be increased by a factor 1.9. 

5. Discussion and conclusions 

We have shown that a gravitational angular-momentum flux, or torque, can be 
measured directly from the mass distribution in spiral disks. Using multiband photometry, 
we have carried out this measurement for Ml 00. The estimated torque implies a dynamical 
evolution time 5 — 10 X 10® yr, with the main source of error being the stellar mass-to-light 
ratio. Our measurement is somewhat sensitive to the orientation of the galaxy, which 
we have estimated from external kinematic data. Dust and young stars appear not to 
be problematic. We obtain consistent results in all three photometric bands. Although 
dust lanes and local enhancements in the surface brightness due to young stars are more 
prominent in the shorter- wavelength images, these features do not affect the m = 2 Fourier 
components or the torque very much, presumably because they cover rather small areas. 

MlOO may not be a typical spiral. It is bright, and its grand-design structure may have 
been stimulated by gravitational perturbations from its small companions. The "bulge" is 
not well flt by a de Vaucouleurs proflle, and it may harbor a small bar [Shaw et al. (1995)]. 
Therefore, the implications discussed in the next paragraph need to be conflrmed by 
applying these methods to a fair sample of spirals. 

We cannot resist the temptation to speculate. If the spiral structure we see today 
is typical of MlOO's past, it has very substantially altered the radial distribution of the 
oldest stars in this galaxy. Whereas the self-similar infall model of Gunn (1981) predicts 
that the disk scale length should increase with time as a consequence of the arrival of 
higher-angular-momentum gas, the torques discussed here would decrease the scale length 
by removing angular momentum from most of the pre-existing stars. It is unknown which 
of these two effects dominates, but the implications for redshift-angular-diameter tests 
based on spirals are not encouraging. More positively, we can say that our results are 
consistent with the speculation of Toomre (1981) that spirals stabilize themselves against 
nonaxisymmetric modes by redistributing mass. Perhaps spirals are enormous accretion 
disks, and because they are so well-resolved, they may offer lessons that can be applied to 
smaller self-gravitating disks such as those believed to exist around protostars. 

We would like to thank James Gunn and James Rhoads for discussions. We are 
indebted also to Edward Fitzpatrick for help with the IRAF image processing software. This 
work was supported by NASA's Astrophysics Theory Program under grant NAGW-2419. 
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Fig. 1. — (a), (b), (c): Images of MlOO in i, r, and g bands, respectively. Scale bar on right 
= 1'. (d) Rectified g image, using i = 27°, PA = 158°. Straight line across image drawn at 
158°. In all panels north is up and east is left. 
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Fig. 2. — Comparison of original data to those reconstructed from the first 10 Fourier 
components, (a) Rectified i image, (b) Fourier reconstruction for i. (c), (d): Like (a) and 
(b), but for g. Scale as in Fig. 1, but note different orientation. 




Fig. 3. — Axisymmetric surface brightness profiles converted to Johnson photometric bands. 
Solid: /. Dashed: R. Dotted: V. Straight line is an exponential with scale length 64". 
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Fig. 4. — Fourier components m = 1, 2, 3, 4 of surface brightness. Left panels are moduli 
normalized to axisymmetric component: (cr^„ + ^smY^'^ I ^o- Right panels are the phases, 
tan~-^(— o"sm/o"cm)- Solid lines: /; dashed: R; dotted: V. 




Fig. 5. — Arm - interarm density contrast. Brightness in magnitudes relative to azimuthal 
mean, plotted against azimuth at 5 radii, after subtraction of odd Fourier harmonics. Solid 
lines: /; dashed: R; dotted: V. 
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Fig. 6. — Rotation curves for MlOO calculated from photometry. Solid line: /; dashed: R; 
dotted: V . Symbols: Velocity data from the the indicated sources; Rubin et a/.data rescaled 
from 2 = 35° to 2 = 27°. 
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Fig. 7. — Torque versus radius for all three bands. Upper panel includes contributions from 
1 < m < 10, lower from m = 2 only. Solid lines: /; dashed: R; dotted: V. 
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Fig. 8. — Sensitivity of the torque to rectification parameters. All curves rescaled by the 
same factor such that correct (solid) curve has peak value unity. Upper panel: artificial 
galaxy. Lower panel: MIOO. For full explanation, see text. 
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Fig. 9. — Evolution time Eq. (15). Solid lines are for /, dashes are for R, and dots are for 
V. 



